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Abstract. The surface critical behaviour of the interacting self-avoiding trail is 
QJ ' examined using transfer matrix methods coupled with finite-size scaling. Particular 

H ■ attention is paid to the critical exponents at the ordinary and special points along the 

collapse transition line. The phase diagram is also presented. 
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1. Introduction 

Lattice self-avoiding walks have been used as models for real, linear polymers in solution 
for over three decades[l]. The quality of the solvent may be introduced by the inclusion 
of short-ranged interactions in the model; typically an attractive energy is included for 
non-consecutive nearest-neighbour occupied lattice sites. This model is the standard 
Interacting Self- Avoiding Walk model (ISAW) or 0-point model[2l[3]. The model has 
been shown to accurately predict the critical behaviour of a wide range of real linear 
polymers in solution, not only in the high-temperature phase, but also at the collapse 
transition, which occurs as the temperature is lowered, at the G temperature. The 
model is successful because it captures the strong entropic repulsion between different 
portions of the polymer chain (the self-avoidance), as well as the effect of the difference of 
affinity between monomer-monomer contacts and monomer-solvent contacts (attractive 
interaction). 

Whilst the relevant physical dimension in polymer physics would usually he d = 3, 
the ISAW model has been much studied in two dimensions. This is partly motivated by 
the realisation that d = 3 is the upper critical dimension of the collapse transition, and 
that the model in two dimensions provides an interesting playground. In this paper we 
shall concentrate on the two-dimensional square lattice. 

One could ask whether the ISAW is a special model, or whether other models which 
include the same basic ingredients would have the same critical behaviour. Two related 
models were introduced to examine this question: the O(n=0) model introduced by 
Nienhuis, which we will refer to in this paper as the Vertex-Interacting Self-Avoiding 
Walk model (VISAW)[1], and the Interacting Self- A voiding Trails model (ISAT)|5]. In 
both of these models the self-avoidance constraint is relaxed in that lattice sites may be 
visited more than once, but the lattice bonds may only be visited once. In the VISAW 
model the walk is not allowed to cross, but in the ISAT model the walk is allowed to 
cross. 

Whilst the VISAW and ISAT models have the same critical behaviour as the 
ISAW in the high-temperature phasejH [6], the situation is different at the collapse 
transition. For the VISAW model, a mapping to an integrable 19 vertex model 
allows the exact calculation of the correlation exponent u = 12/23 and the exponent 
7 = 53/46 [7], as compared to z/ = 4/7 and 7 = 8/7 for the ISAW at the collapse 
transition [8]. The situation for the ISAT model is far less clear. Early studies found 
a number of contradictory results [9l [TOl [11], [12]. Some authors claimed that the 
results were compatible with the ISAT model and the ISAW model being in the same 
Universality class at the collapse transition [10], whilst others found results which were 
incompatible im [T2] . Recently this model was re-examined using transfer matrices |13]. 
and the conclusion arrived at was that the transition was similar to the collapse 
transition in the VISAW model (they have the same correlation length exponent u but 
different 7 exponents). The collapse transition for both the VISAW and ISAT models 
is clearly not in the usual 0-point transition; the collapse transition in these models has 
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an extra transition line leaving it separating two distinct dense phases. 

In our previous study of the ISAT model|13], we gave results for the bulk-critical 
exponents and showed that these results were different from the results of Owczarek and 
PrellbergjTT]. We argued that this was due to a breakdown of the usual identification 
of the u exponent with the geometrical exponent derived from the radius of gyration. 
In the current paper, we extend our study by introducing attractive interactions with 
a surface, and examine the surface critical behaviour of the model. In their paper [TT]. 
Owczarek and Prellberg also presented some results for the surface critical behaviour. 
Some of their results corresponded to explicit calculations using Monte Carlo, and others 
were derived by plausibility arguments. Our results concord with their explicit results, 
but disagree with the others. The surface critical exponents we find are consistent with 
the bulk results found in [T3l. 



2. Model and Transfer Matrix Calculation 

The ISAT model studied here is defined as follows: consider all random walks on the 
square lattice which do not visit any lattice bond more than once. Doubly visited 
sites may correspond to either crossings or "collisions" ; both are assigned an attractive 
energy —e. The walk is allowed to touch, but not cross, a surface defined as a horizontal 
line on the lattice. Each step along the surface is assigned an attractive energy —es- 
For the transfer matrix calculation that follows it is convenient to consider a strip of 
width L with an attractive surface both sides of the strip. This is not expected to 
change the behaviour in the thermodynamic limit L — )■ oo; the bulk critical behaviour 
should not depend on the boundary conditions and when calculating the surface critical 
behaviour, a walk adsorbed to one surface needs an infinite excursion in order to "see" 
the other surface. Additionally, the finite-size scaling results which link the eigenvalues 
to of the transfer matrix to the scaling dimensions x* and x^ (see flTSj) and f lT6|) ) rely on 
the conformal mapping of the half plane (one adsorbing surface) onto a strip with two 
adsorbing surfaces[14j. A typical configuration for the ISAT is shown in Figure [H 
The partition function for the model is 

Z=Y^ K'^u^H'"', (1) 

walks 

where K is the fugacity, Ug = exp{l3es) and r = exp(/3e). A^ is the length of the walk, 
Ns is the number of steps on the surface, and Nj is the number of doubly-visited sites. 
The average length of the trail is controlled by the fugacity K through 

(-)^A-^. (2) 

As K increases from zero, (A^) increases, diverging at some value K = K*{us,t). To 
start we consider what happens in the absence of the adsorbing boundary. For r small 
enough, 

{N)^{K*{oo,,t)-K)-\ (3) 
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Figure 1. A self-avoiding trail model showing the vertex crossings and the vertex 
collisions, both weighted with a Boltzmann factor r. Surface contacts are weighted uj^ 
and a fugacity K is introduced per walk step. The trail is shown on a strip of width 
L = 5. 



K 



(Oc 



whilst for large enough r the divergence is discontinuous. Whilst (A^) is finite, the 
density of occupied bonds on an infinite lattice is zero, whilst once (A^) has diverged 
the density is in general finite. For small enough r the density becomes non-zero 
continuously at i^ = K* and for large enough r the density jumps a.t K = K*. K* may 
then be understood as the location of a phase transition, critical for r < Tcou and first 
order for r > Tco\\. The problem of infinite walks on the lattice is equivalent to setting 
K = K* and varying r, then it may be seen that for r < Tcou the density is zero and is 
non-zero for r > Tcow- It then follows that Tcoii defines the collapse transition point. 

Now let us consider the effect of the adsorbing boundary at constant r. For Us 
small, the entropic repulsion of the wall is strong enough for the walk to remain in the 
bulk. Once Ug is large enough for the energy gain to overcome the entropic repulsion, 
the walk will visit the boundary a macroscopic number of times, and the walk adsorbs to 
the surface. These two behaviours are separated by a;^ = w*. For Us < ujg the behaviour 
of the walk is not infiuenced by the wall, and K* is independent of Ug- The transition 
K = K* if critical (r < Tcou) corresponds to ordinary critical behaviour. However, for 
Ws > w*, K* is a function of w^, and the transition is referred to as a surface transition. 
The point K = K*, Ug = ou* is referred to as the special critical point (again r < Tcou). 

As the critical value K* is approached, and in the absence of a surface, the partition 
function ([T]) and the correlation length ^ are expected to diverge, defining the standard 
exponents 7 and u: 

^r^lK-K*]-" (4) 

Zr^\K-K*\-^ (5) 

The effect of the surface on the walk is to introduce an entropic repulsion, pushing the 
walk away from the surface. The number of allowed walks is reduced exponentially if 
the walk is constrained to remain near the surface, in particular if one or both ends of 
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the walk are obliged to remain in contact with the surface. In this case the divergence 
of Z is modified, and two new exponents are introduced, 71 and 711. Defining Zi and 
^11 as the partition functions for a walk with one end, and both ends, attached to the 
surface respectively, then: 

Zi r^ \K - K*\-^' (6) 

Zn -- \K - K*\-'^'' (7) 

Whilst the bulk exponents, such as u and 7, are the same at an ordinary critical point 
and at the special critical point, the surface exponents 71 and 711 differ. The exponents 
z/, 7, 7i and 711 are related by the Barber relation [TSj: 

I/ + 7 = 271 -711. (8) 

This partition function may be calculated exactly on a strip of length L^ — )■ 00 and 
of finite width L by defining a transfer matrix T. If periodic boundary conditions are 
assumed in the x-direction, the partition function for the strip is given by: 

Zl= lim Tr(r^-). (9) 

The free energy per lattice site, the density, surface density and correlation length for 
the infinite strip may be calculated from the eigenvalues of the transfer matrix: 

fiK,u,,T) = ^\ni\o), (10) 

p(^,-s,r) = -^^, (11) 



LAo dK 

's_<- 

Ao du.. 



ps{K,u.,r) = ^p, (12) 



e(ir,a;„r)= In 



(13) 



, aJ ' 

where Aq and Ai are the largest and second largest (in modulus) eigenvalues. 

Our first task is to find estimates of K*{uJs,t). For this, two distinct methods are 
used, which we now describe. 

(i) For K < K* the largest eigenvalue of the transfer matrix T corresponds to the 
empty lattice; the length of the walk is finite, but the lattice strip is infinitely 
long, and so the probability of finding a non-empty column is zero. The largest 
eigenvalue is then Aq = 1. The divergence of the walk length is identified with the 
value K = Kl for which ^ — )• 00. This occurs when Ai = Aq = l[16j. 

(ii) An estimate for the critical point where the length of the walk diverges may be found 
using phenomenological renormalisation group for a pair of lattice widths [T7j, L and 
L'. The estimated value of K* is given by the solution of the equation: 

^ = ^ (14) 
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Both these methods give finite-size estimates KI{us,t) which should converge to the 
same value in the limit L — )■ oo. 

The critical dimensions of the surface magnetic and energy fields may be calculated 
from the first few eigenvalues of the transfer matrix: 

Lin ^ 

(15) 



x„ 



TC 



Lin 

x" 



vr 



(16) 



with A2 the eigenvalue with the third largest absolute value. 

The surface scaling dimensions x^ and x| may be related to the surface correlation 
length exponent Ug and the exponent r^y, controlling the decay of the correlation function 
along the surface, through standard relations 

^s = ^^, (17) 

1 — x| 

^11 = 2<. (18) 

The entropic exponent 711 is related to ri\\ through: 

7n = Kl-^||)- (19) 

For a more detailed discussion of the transfer matrix method, and in particular how 
to decompose the matrix, the reader is referred to the article of Blote and Nienhuis [1]. 

3. Results 

For clarity, we will present separately the results found setting Ai = 1 and those found 
using the phenomenological renormalisation group equation ( IT^ . 

The finite-size results obtained are, where possible, extrapolated on the one hand 
using the Burlisch and Stoer (BST) extrapolation procedure fl8| and on the other hand 
fitting to a three point extrapolation scheme, fitting the following expression for quantity 

Xl = X^ + aL-\ (20) 

Calculating X^o, a and b require three lattice widths. The extrapolated values Xoo 
clearly will still depend weakly on L, and the procedure may be repeated, however weak 
parity effects can be seen in their values, impeding further reasonable extrapolation by 
this method. 

3.1. Results setting Ai = 1 

In Figure [2] we present the curves of Kl{us, r) as a function of cj^ for r = 0. The case 
where r = corresponds to the pure self-avoiding walk in the presence of an attractive 
surface, already studied using transfer matrices p^ 120] . The interacting self-avoiding 
trail maps onto the kinetic self-avoiding trail model when r = 3. It was conjectured 
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Figure 2. KI{lUs,t) calculated setting Ai = 1 for r = 0. This case corresponds to 
the adsorption of the standard SAW model. 
0.5 




CO 



Figure 3. K1{uJs,t) calculated setting Ai = 1 for r = 3. r = 3 is expected to 
correspond to Tcou, corresponding to the collapse transition line for Us < w*. 
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that the kinetic self-avoiding trail corresponded to the collapsing self-avoiding trail, 
leading to the identification Tcoh = 3[llJ. This identification seems to be verified by 
the results presented in references |ll] and [13]. The curves of KI{ll!s,t) for r = 3 
are shown in Figure |3l The values calculated from the intersections of successive 
curves are shown in Table [1] which places the special surface transition for r = 3 at 
K* = 0.3333 ± 0.0001, w* = 2.427 ± 0.002. We calculated the finite-size estimates 
for the density pl{K*,u)*). These estimates can be seen to have a finite limit of 
Poo = 0.35 ± 0.02. In Table [Hwe also give this reciprocal fractal dimension calculated 
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Figure 4. Phase diagram in the K = K1{lUs,t) plane found setting A = 1 and 
identifying the adsorption transition with the the solutions oi Kf^{iOs, r) = K1_^-^(ujs, r). 
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at the special transition for r = 3. Whilst the question of the density at the special 
point in the infinite lattice limit remains open, given the lattice widths used in our 
calculation, we can clearly see that we recover a geometrical exponent compatible with 
the results Owczarek and Prellberg[Tl]. This indicates that the thermal exponent i' 
and geometrical exponent 1/df are likely to be different for this model at the collapse 
transition. 

For both r = and r = 3 the curves may clearly be seen to cross at a point, 
which defines the adsorption transition ^^(r). By varying r we may map out the 
adsorption line, and hence find the phase diagram in the K = K*{ujs, t) plane. This is 
shown in figure HI The vertical line, corresponding to the collapse transition, has been 
added by hand at Tcoh = 3; the boundary is not expected to influence the location of 
the collapse transition until the walk adsorbs to the boundary; the surface interaction 
below u* will have a finite perpendicular correlation length ^^ associated with it. In the 
thermodynamic limit, the bulk of the lattice will not be influenced surface, and bulk 
collapse should be unaffected. It is only at u*, were ^j^ — )■ oo, that the bulk will be 
influenced by the surface, and the walk will adsorb to the surface. This picture has 
been proven for the case of the adsorption of collapsing self-avoiding polygons in three 
dimensions [21] . 

Having set r = 3 using the conjecture from the kinetic SAT model, there are 
different methods for estimating the location of the special surface transition, where 
collapse and adsorption occur simultaneously. The simplest is to use the crossings 
of K^ to estimate the adsorption point w*, results for which are shown in Table [H 
However, we may also use the scaling behaviour of the order parameter to set up 
a phenomenological renormalisation scheme. An order parameter for the adsorption 
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Table 1. Location of the special point calculated by setting Ai = 1 and solving 
K'^{u)*,t) = K'^^-^^{ui*, t) for T = 3. The results have then been extrapolated using the 
BST extrapolation scheme. The reciprocal fractal dimension l/d/ is given, calculated 
from the densities pL and pl+i- 



K* 



PL+l 



Udf 



3 


2.220342 


0.340394 


0.501949 


0.565358 


4 


2.267276 


0.337428 


0.474241 


0.545825 


5 


2.296043 


0.336003 


0.455622 


0.534034 


6 


2.315608 


0.335210 


0.442176 


0.526218 


7 


2.329838 


0.334723 


0.431958 


0.520713 


8 


2.340687 


0.334403 


0.423892 


0.516664 


9 


2.349254 


0.334182 


0.417336 


0.513588 


10 


2.356203 


0.334022 


0.411881 


0.511191 


11 


2.361962 


0.333904 


0.407254 


0.509286 



BST cx) 2.4267 



0.3333 



0.3539 



0.495 



Three point extrapolated results 



3 


2.4227 


— 


— 


— 


4 


2.4250 


0.333267 


0.327394 


0.488281 


5 


2.4262 


0.333291 


0.332348 


0.490432 


6 


2.4271 


0.333304 


0.334328 


0.491926 


7 


2.4277 


0.333312 


0.334880 


0.493045 


8 


2.4281 


0.333317 


0.334636 


0.493983 


9 


2.4284 


0.333321 


0.333949 


0.494649 



Table 2. Location of the special point calculated by setting Ai = 1 and solving 
^L.L+ii^t^''') — YL+i,L+2{^t,T) = y^ for T — 3. The results have then been 
extrapolated using the BST extrapolation scheme. The precision of the estimates 
for 3^oo is not sufficient for the BST extrapolation to be very accurate. The BST value 
is yoc ~ 1-08 which is a little low. Plotting the points, a reasonable estimate would be 
y^o = 1.12 ±0.05 



L 


^*s 


K* 


yL = {i- <i>s)lv 


3 


2.644273 


0.306760 


0.635164 


4 


2.612544 


0.315671 


0.675046 


5 


2.587753 


0.320837 


0.714237 


6 


2.568090 


0.324099 


0.751030 


7 


2.552274 


0.326280 


0.784704 


8 


2.539271 


0.327805 


0.815357 


9 


2.528203 


0.328916 


0.843678 


10 


2.519807 


0.329713 


0.866822 


BST oo 


2.431 


0.3339 


see text and caption 
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Figure 5. Plot of the order parameter O = {Ns)/{N) for (a) r = 3 and (b) t — 4, 
calculated setting A = 1. 

(a) 




(b) 
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transition may be defined as 



(Ns) 
(N) ' 



(21) 



For a continuous adsorption transition, the order parameter O vanishes for Ug < cu* 
whilst for Us > ujg O is finite. At the transition point the order parameter is expected 
to scale as : 



Ol{uJs) = L-^'-^^y-'d [ly^ius - oul)] . 



(22) 
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Defining Yl^l, by 

\n[0,^{u,)/0,ius)] .... 

^^'^'^^^) = h[[Lji7] ' ^^^^ 

we may set up a phenomenological RG scheme for estimating tlie adsorption transition, 
as well as the exponent (pg. These are given by the solutions of the equation 

>L,L+lK) = >L+1,L+2K) = yi- (24) 

where y^ is the finite size estimate of the exponent 3^oo = (1 ^ (l>s)/^- 

Plots of the order parameter are shown in Figure [5] and results for yi are shown 
in Table [2l The values y^ are far from their asymptotic values, and the extrapolation 
using BST is not very conclusive. The best estimate from BST is yi = 1.08. Plotting 
this value along with the data points shows it to be a low estimate. A reasonable 
extrapolation of the points seems to be y^ = 1.12 ± 0.05, although these values should 
be taken with care. 

3.2. Results from Phenomenological RG 

In the previous section we calculated our estimate K*{us,t) by setting Ai = 1. This 
gives an overestimate of K*. Calculating K* using equation [H] gives more asymptotic 
results, but at the cost of needing two lattice widths, reducing the number of available 
data points for extrapolation. Additionally, the phenomenological RG method used 
identifies both the ordinary and the special surface transitions. The phase diagram in 
the {K, Us) plane is shown in figure [6] for r = 3. It may be seen that there are two 
points on the phase diagram where the lines seem to cross at (or close to) a point. 
These two points are the fixed points corresponding to the ordinary and the special 
surface transitions. The location of these points may be estimated using a three width 
phenomenological RG scheme. From equations (fT3|l . (fT4|) . (fTSll and (fT8|) it may be 
seen that we may identify the fixed points with crossings of r]\\ , and find corresponding 
estimates for r^y . The plots of r^y as a function of Ug for r = 3 are shown in Figure U\ 
whilst the phase diagram in the r, Ug plane whith K = K* is shown in Figure [HI The 
estimates for the location of the fixed points for the ordinary and special transitions, 
along with the corresponding exponent estimates, are shown in Table [3] and Table HI 

The results presented in Table [3] and Table H] are calculated setting r = 3. The 
identification t^ow may be checked by looking for the special transition fixed point in the 
full phase diagram at the cost of using an extra (fourth) lattice width. Unconstrained 
results are shown in Table [5l Extrapolated values of r give Tcoii = 3.000 ± 0.001, which 
is consistent with the identification used elsewhere. 

In the results presented above we have used ( IT^ . or equivalently (TT5l) . to identify 
the fixed points. However the correlation length defined using the largest and third 
largest eigenvalues also diverges in the thermodynamic limit, and so we may use a 
phenomenological RG based on Xe using ( TT6|) . This leads to estimates for Vg-, the 
correlation exponent along the surface as well as alternative estimates for the location 
of the fixed points. The results are shown in Table El 
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Table 3. Ordinary point location for t = 3 and estimates for ?7Jj"^'^. 
extrapolations are shown in the second half of the table. 



The three point 



K* 



m 



3 


0.494271 


0.336986 


1.734977 


4 


0.369013 


0.334303 


1.852456 


5 


0.313660 


0.333636 


1.892394 


6 


0.289337 


0.333451 


1.907045 


7 


0.277118 


0.333387 


1.913419 


8 


0.270290 


0.333362 


1.916570 


9 


0.266155 


0.333351 


1.918282 


10 


0.263473 


0.333339 


1.919288 



BST cx) 0.256729 0.33333 1.92133 



Three point extrapolated results 



3 


0.225177 


0.333230 


1.931118 


4 


0.256335 


0.333488 


1.921191 


5 


0.256947 


0.333401 


1.921116 


6 


0.256609 


0.333339 


1.921221 


7 


0.256325 


0.333339 


1.921278 


8 


0.255885 


0.333409 


1.921363 



Table 4. Special point location for r = 3 and estimates for r;? 
extrapolations are shown in the second half of the table. 



The three point 



L 



K* 



sp 



3 


2.462062 


0.332813 


-0.110939 


4 


2.451809 


0.333074 


-0.105536 


5 


2.446565 


0.333181 


-0.102198 


6 


2.443407 


0.333234 


-0.099858 


7 


2.441324 


0.333264 


-0.098106 


8 


2.439855 


0.333283 


-0.096728 


9 


2.438768 


0.333294 


-0.095606 


10 


2.437932 


0.333302 


-0.094666 



BST oo 2.43245 0.33331 



-0.08449 



Three point extrapolated results 



3 


2.434491 


0.333327 


-0.086941 


4 


2.433346 


0.333328 


-0.083971 


5 


2.433139 


0.333332 


-0.083208 


6 


2.432852 


0.333340 


-0.082012 


7 


2.432683 


0.333317 


-0.081190 


8 


2.432473 


0.333339 


-0.079932 
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Table 5. Special point location with r unconstrained, using four lattice widths. 
Estimates for rjl^ are also shown. The three point extrapolations are shown in the 
second half of the table. 



L 


T* 


^t 


K* 




3 


3.132743 


2.514601 


0.329697 


-0.123477 


4 


3.074773 


2.486484 


0.331186 


-0.115417 


5 


3.048363 


2.472046 


0.331897 


-0.110607 


6 


3.033676 


2.463154 


0.332307 


-0.107244 


7 


3.024799 


2.457278 


0.332562 


-0.104764 


8 


3.019029 


2.453143 


0.332731 


-0.102840 


9 


3.015082 


2.450106 


0.332848 


-0.101298 


BST 


cx) 2.99967 


2.440 


0.333291 


-0.0874 




Three point extrapolated : 


results 




3 


3.003123 


2.438430 


0.333265 


-0.091790 


4 


2.997885 


2.432266 


0.333423 


-0.084854 


5 


3.000002 


2.433993 


0.333345 


-0.086107 


6 


3.000077 


2.433675 


0.333338 


-0.085025 


7 


3.000388 


2.433813 


0.333325 


-0.084784 



Table 6. Special point location with t = 3 by applying the phenomenological RG 
method using x^., leading to estimates for Ug = 1/(1— x^). The three point extrapolations 
are shown in the second half of the table. 



L 



K* 



3 


2.380861 


0.334825 


1.352725 


4 


2.423360 


0.333648 


1.309484 


5 


2.439110 


0.333307 


1.291723 


6 


2.445939 


0.333198 


1.283136 


7 


2.449061 


0.333167 


1.278793 


8 


2.450430 


0.333166 


1.276705 


9 


2.450901 


0.333175 


1.275924 


10 


2.450890 


0.333188 


1.275943 



BST oo 2.4512 



0.3332 



1.275 



Three point extrapolated results 



3 


2.456846 — 


1.267398 


4 


2.454952 — 


1.268746 


5 


2.453232 — 


1.271502 


6 


2.451999 — 


1.273811 


7 


2.451239 — 


1.275283 


8 


— — 


1.276304 
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Figure 6. KI{uJs,t) calculated using phenomenological RG for t — 3. 
0.4 1 . 1 . 1 . 1 . 1 . 



0.35- 



K 0.3- 



0.25- 




Figure 7. Exponent 7711 from crossings of ^/L for r = 3 and L' = L + I 




In Table [2] we calculated estimates, yi, of 3^oo = (1 ~ 4's)/^ from (l23l) and ^ 
calculating K*{ujs, r) by setting Aq = 1. In Table [7] we present the analogous calculation 
but with K* calculated using phenomenological RG. 



3. 3. Summary of results 

The transfer matrix calculation gives the exponent r]\\ directly. At the ordinary transition 
we find a value ?7|j"''^ = 1.9213 ± 0.0001. Using the results conjectured in [13], z^ = 12/23 
and 7 = 22/23 with ^ and ^, we find 7°'''^ = 0.499 and -f°f = -0.481. These values 
agree with the value 7™'^ calculated with Monte Carlo in yjj, but not with the value of 
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Figure 8. Phase diagram calculated by the crossings of t^h . 
4i 1 1 . 1 . 1 1 1 



3.5- 
3- 



«2.5h 



2F 

1.5 



Adsorbed Phase 




p=0 phase 



Collapsed Phase ■ 







Table 7. Location of the special point calculated by fixing K* by phenomenological 
RG and solving Yl{lu*,t) = Yl+i{uj*,t) for r = 3. The results are have then been 
extrapolated using the EST extrapolation scheme. 



L 


K* 


^: 


Yl 


3 


2.675620 


0.326891672428 


0.715973 


4 


2.622565 


0.329233649135 


0.767296 


5 


2.587417 


0.330544853438 


0.809961 


6 


2.562356 


0.331342978422 


0.845875 


7 


2.543648 


0.331858273706 


0.876360 


8 


2.529207 


0.332206642617 


0.902462 


9 


2.517800 


0.332450587489 


0.924981 


10 


2.507957 


0.332633569013 


0.945751 


EST oo (3-8) 0.333377 


2.427 


1.19 


Three point extrapolated results 


3 


0.330183 


2.348631 


— 


4 


0.331123 


2.369608 


— 


5 


0.333704 


2.388620 


— 


6 


0.333557 


2.401071 


— 


7 


0.333463 


2.413110 


— 


8 


0.333647 


2.222366 


— 
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7°5^*^. It should be noted that the values were obtained using u = 12/23 . At the special 
transition we find r]f^ = —0.085 ±0.003, leading to 7^^ = 1.022 and 71]' = 0.566. Again, 
7i^ is consistent with the value given in [H], but not 7^^, and again we used the value 
z/ = 12/23. In both cases in [11] the values of 711 are calculated using the assumed 
values of 7 = 1 and z/ = 1/2, which we believe not to be the correct values. 
We calculated 0^ in three ways: 

(i) From the result z/, = 1.275 ±0.002 and using u = 12/23 leads to 0^ = u/us = 0.414. 
(ii) Using 3;^ = (1 - 0,)/zy = 1.12 ± 0.05 from Table El giving 0, = 0.405. 
(iii) Using 3;^ = (1 - 0,)/z/ = 1.19 ± 0.05 from Tabled giving 0, = 0.379. 

In all three cases the result is close but a little lower than the value (ps = 0.440 ± 0.010 
given in [H] , but the determination of this exponent is the least accurate of the results 
presented here. 



4. Discussion 

In this paper we have investigated the surface critical behaviour of the Interacting 
Self-Avoiding Trail. The results confirm the identification of the collapse transition 
as occurring at i^cou = 1/3 and Tcou = 3 and locates the adsorption transition for 
''"coll = 3 at cjs = 2.45 ± 0.05. For the standard Interacting Self-Avoiding Walk model 
an additional line has been observed in the collapsed phase, separating a region where 
the collapsed walk is adsorbed to the surface from a region where the walk is desorbed 
from the surface [22]. This transition corresponds to a wetting of the surface by the 
perimeter of the collapsed globule, and so does not correspond to a singularity in the 
bulk free energy. Whilst the method used in this paper does not naturally throw up 
this transition line, an a-posteriori reexamination of the order parameter defined in ( 12T|) 
for the standard ISAW model does show a signature of this transition (see the curves in 
ref [20]). Such a signature seems to be absent in for this model, see Figure |5] for plots 
with r = 4 > Tcoii- Whilst this should not be taken as a strong indication, it would 
be of interest to examine whether such a line does exist in this model, or whether the 
particular type of collapse suppresses this transition. 

In conclusion, in the article we use finite-size scaling and transfer matrix methods 
to investigate the surface critical behaviour of the interacting self-avoiding trail model. 
Using the conjectured results for u and 7 from reference |13] (z/ = 12/23 and 7 = 22/23) 
we find values of 7°'''^ and 75*^ in agreement with the numerical results of Owczarek and 
Prellbergjll]. We do not, however, find agreement with the conjectured values of 7°^^^ 
or 7^2*. The values of v and 7 needed for this agreement are different from the values 
u = 1/2 and 7 = 1 proposed in [llj. 
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